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SUMMARY 


The formulas of supersonic wing theory for source, doublet, and 
vortex distributions are reviewed and a systematic presentation is 
provided which relates these distributions to the pressure and to the 
vertical indu ced velocity in the plane of the wing. It is shown that 
care must be used in treat ing the singular it ie s involved in the analysis 
and that the order of integration is not always reversible. Further, it 
is shown that the use of the complex variable can often facilitate the 
calculation of the integrals involved. Certain special applications are 
included to illustrate the concepts presented. 


INTRODUCTION 


One of the most fundamental approaches to the analytical investiga- 
tion of linearized wing theory, throughout the subsonic and supersonic 
Mach number range, stems from the use of certain elementary mathematical 
expressions which are identified physically with sources, doublets, and 
vortices in the fluid medium. By means of these expressions, boundary- 
value problems involving wings with thickness, camber, and angle of 
attack can be solved. These problems are divided into two categories: 
one, involving symmetrical bodies with thickness and no lift, is analyzed 
by means of source distributions; and the other, involving lifting plates 
without thickness, is analyzed by means of doublet and vortex distribu- 
tions. 

All these distributions require the treatment of singularities in 
the mathematical analysis. Thus, for subsonic Mach numbers, the concept 
of Cauchy^ principal part plays an important role in the calculation of 
integrals arising in the development of lifting— line and two-dimensional 
section problems. In supersonic wing theory, the Cauchy principal part 
is again used in the treatment, for example, of conical— flow problems as 
in reference 1, but, because of the Mach lines and cones appearing in 
the physical flow and in the hyperbolic geometry of the differential 
equation, other techniques in handling improper integrals are needed. 



2 


NACA TN 2252 


The integrals in supersonic wing theory thus require, in general, more 
careful attention to the discontinuities in the integrand and, as an 
illustration, indiscriminate use of such standard devices as inversion 
of the order of integration may lead to incorreet results. 

When problems of the first kind are involved, that is, when pre- 
scribed distributions are to be integrated (as for the problem of finding 
the pressure on a wing with symmetrical thickness), a guide to the proper 
method of calculation is often furnished by physical intuition. However, 
when problems of the second kind arise, that is, problems, the solution of 
which depends upon the inversion of an integral equation (as the flat 
plate of arbitrary plan form), the mathematical methods are more abstract. 

The purpose of the present report is: first, to review the formulas 
of linearized wing theory in which source, doublet, and elementary- 
vortex distributions are introduced and to relate these distributions to 
the pressure and to the vertical induced velocity in the plane of the 
wing; second, to show that the use of the complex variable can often 
facilitate the calculation of the integrals involved; and finally, to 
present certain special applications which will illustrate the basic 
concepts. . 


LIST OF IMPORTANT SYMBOLS 


Mo 

<1 


q 

T r 




pressure coefficient 
free— stream Mach number 

loading coefficient (pressure on the lower surface minus pressure 
on the upper surface divided by free— stream dynamic pressure) 


free— stream dynamic pressure 


(2 p ° T ° 2 ) 


hyperbolic distance between points x,y, z and xi,yi,0; 
r c =J (x-xi) 2 -f3 2 (y-yi) 2 -§ 2 z 2 


V 0 velocity of the free stream 
x,y, z Cartesian coordinates 
r,s,z characteristic coordinates 
5,q,z oblique coordinates 

mi cotangent of the angle between the q and x axes 
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m 2 cotangent of the angle "between the | and x axes 

M-i a/ 1+mi 2 

H2 y l+ms 2 

u perturbation velocity in x direction 

w perturbation velocity in z direction 

p y^-i * 

X slope of stream surface (w/V Q ) 

cp perturbation velocity potential 

A Jump in value of the quantity considered across the z = 0 plane 

Sub script 

u value of a quantity on the upper surface of a wing (z = 0 plane) 


GENERAL THEORY 


It Is well known that the analysis of thin wings at supersonic 
speeds and at small singles of axtack can be expressed in mathematical 
terms as a boundary— value problem for the wave equation. If ft repre- 
sents a velocity potential or any one of the velocity components them- 
selves, this equation can be written 


P 2 — ft yy — ^zz = ® (l) 

where the z = 0 plane is the plane of the wing, the free— stream veloc- 
ity V 0 is directed along the x axis and P 2 = Mq 2 — 1, being tbe 
free— stream Mach n umb er. 

Of the many ways of solving the boundary— value problems associated 
with the wave equation, the most convenient for the present purpose is 
the Volterra solution. In reference 1, a discussion was given of the 
application to aerodynamic problems of Volterra* s method. Thus, the 

solution for ft in terms of A ft and A^> the Jump in ft and its 

gradient in crossing the z = 0 plane can be written, if 
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r c -J (z-xi) 2 -p 2 (y-yi) 2 — (3 2 z 2 


( 2 ) 


in the form 


ft 


•-h “■= ° osi1 


dft 


z— z n 


P//(y-7i) 2 + z 2 


dzi dyi + 


JL_ _d_ 

,2n ^x 



^ Aft 


z(z-zi) dzi dyi 
E(y-yi ) 2 + z 2 ] r c 


(3) 


where the region t is that portion of the z = 0 plane ‘lying within 
the forecone from the point x,y,z (the forebranch of the hyperbola 
r c 2 .=0). 

Equation (3) represents a general solution to the wave equation and 
has yet to be put in a form which represents directly a solution to a 
problem arising in the study of wings. When equation ( 3) is so adapted 
in the nezt section, it will represent the velocity potential due to a 
distribution of sources, doublets, and elementary horseshoe vortices, 
the strength of which are given in terms of the wing shape and loading. 

The following study of the adaptation of equation (3) to the partic- 
ular boundary— value problems of wing theory requires the introduction of 

two axial systems other than 
the z,y,z Cartesian coordi- 
nates already defined. First, 
the coordinate system 

is defined so that z is normal 
to the plane of the wing while 
I, tj are both normal to z 
(i.e., lie in the plane of the 
wing) and make arbitrary angles 
with the x,y system (see sketch) 
If 

PX = a / l+mx 2 , Ps = y ”i+m 2 2 (k) 

The equations which relate the 
|,tj,z to the z,y,z system 
are 
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Second, the r,s,z coordinate 
system Is also defined so that 
z Is normal to the plane of 
the wing while r,s are both 
normal to z and lie along 
the traces of the Mach cone 
emanating from the origin 
( see sketch) . It is apparent 
that the r,s,z coordinates 
are a special case .of the |,T],z 
system formed when mi=me=P . 
(Notice also that the x,y,z 
coordinates are obtained from 
z when mi=0 and me= ». ) 
The equations which relate 
r,s and z to x,y and z 
are 


I 
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x = j£ (s+r) 


7 = “ (b-t) 
Mq 


z = z 


r = ^ (x-Py) 


b - | (*tfy) > 


z = z 


(8) 


and the Jacobian is 


2 P 

dx dy = —5 dr ds 
Mo 


(9) 


The value of r c is transformed by the equation 


r c =y (^ 2 /Mo 2 ) (r— ri) (&-s x ) H3 2 z £ 


(10) 


TPrEi Thkhih! FDHDAMEKTAL FORMULAS 


As has been indicated already, the next purpose is to relate equa- 
tion ( 3) to the three fundamental formulas arising in wing theory which 
are those relating the velocity potential to Bource, doublet, and vortex 
distributions. These distributions can be expressed in terms of the 
discontinuity in either cp or its gradients in the plane of the wing. 

In this way, the expression for a source distribution is obtained when 
the perturbation velocity potential *P is an even function with respect 
to the z = 0 plane and is expressed as a double integral involving 
A( d<p/dz) where the A notation denotes the jump in the value of 
ckp/dz in crossing this plane; that is 

(IsA -(* 2 .) 

bz \8z / z=0+ \dz / 2=0— 

Similarly, a distribution of doublets is obtained when is an odd 
function with respect to the z - 0 plane and is expressed as a- double . 
integral involving A 9. Finally, a vortex distribution results when «P 
is an odd function and given as a double integral involving the loading 
coefficient Ap/q which is, in turn, equal to 2A<P X /V 0 . In these 
formulas A(5cp/dz), Acp, and Ap'/q are, respectively, the strengths 
per unit area of the sources, doublets, and elementary horseshoe vortices. 
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The Source Distribution 


The velocity potential at the point x,y,z of a unit source at the 
point xi,yi,0 is given by the equation cp = — l/2itr c . The manner in 
which a distribution of these sources affects the potential vill new be 
developed from equation (3). If, in equation (3), ft is set equal to 
the velocity potential, and the potential in turn is assumed to be 
symmetrical above and below the x,y plane (as in the case of a symmet- 
rical airfoil at zero lift), then equation (3) becomes 



A ^ arc cosh — ==== Xl 

dz Py (y-yi) 2 +z 2 


(Dti dyi 


Since the inverse hyperbolic term vanishes on the Mach forecone, the 
partial derivative can be carried through the double integral sign and 
there results 


cp (x,y,z) = “ J d yi (^-) 

where w u is the vertical induced velocity on the upper side of the 
z = 0 plane and A(d<p/3z) = 2w u by reasons of symmetry. Equation (ll) 
is the familiar equation for the velocity potential due to a distribution 
of sources in the xy plane. 

In the coordinate system equation (11) becomes 

* c s sis r rhLkhl aSx dlu (M) 

M-s J T J r c 

T i 

where the area t is transferred to the plane, and r c is given 

in these coordinates by equation (7). 


The Doublet Distribution 


The velocity potential at the point x,y, z of a doublet at the 
point xi,yi,0 is given by the equation cp= j3 2 z/2rtr c 3 . The effect on 
the potential of a distribution of these doublets is not so obvious as 
it was in the case of the sources and considerable care must be shown 
in the development. 
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If the vertical induced velocity is considered always to he equal 
above and below the xy plane, but the potential to be discontinuous 
across it, then equation (3) for the velocity potential becomes 


1_ b_ r /^ z(x-xi)Acp (xi,yi)d3Ei dyi 
E(y-yi) 2 +z 2 ] r c 


In this case the integrand does not vanish at the Mach cone and the 
partial derivative cannot be moved directly through the double integral 
sign. Writing in the limits of integration so that the first integration 1 
is made with respect to yi, 


1 8 pxH3z p 

tei X 471 


^2 


z(x-xi) A cp (xi,yi) 


T i [(y-yi) 2 +z a ] J ( x-xi ) 2 — P 2 ( y— y x ) 2 -p 2 z 2 


(13a) 


where 


and 


Ti = y - p/^x-xi ) 2 -P 2 z 2 


f 2 - 7 + j 7(x-xi) 2 -£ 2 z 2 


If in the yi integral xi is replaced by the value x-f3'z, the 
result is indeterminate. Such an indeterminate form can be evaluated, 
however, by excluding the limit x-Pz from the ‘area of integration. 
Hence consider the integral 


<P 


lim 
e — >0 




z(x-xi) Acp (xi,yi) 

[(y-7i) 2 +z 2 ] J ( x— Xi ) 2 -P 2 ( y— y 1 ) 2 -P 2 z 2 


1 Since in what follows the order of integration is important, the nota- 
tion will be adopted that / dy / dx f(x,y) = / C/f(x,y)dx Jdyj that is, 
the integration is made first with respect to x. When the notation 

// f(x,y)dx dy is used, the order of integration is immaterial. 
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By formal, manipulation 


cp 


z P7 +e \/ z2+g2 A cpCz-^yz^e 2 , yi)dyi 
€ 0 23t ^y-e [(y-yi) 2 +z 2 ] y^-Cy-yi) 2 


1 px-g J z 2 +e 2 ^ r Ta a 

° 1 s x, ^ 


'Yl 


z(x-xi) Acp (xuyi) 


[ (y-yi^+z 2 ] 7 ( x— xi ) 2 -P 2 ( y— y 1 ) 2 -f3 2 z £ 


(13b) 


Application of the mean value theorem to the first integral in equa- 
tion (13b) yields 


lim —a/ z a +e 2 A <p(x-P a/z^+e 2 , y+0s) P 

2rt J^ 


dyi 


y-e [ ( y-yx ) 2 +z 2 3 y^My-yi) 2 


which becomes 


.SSb 


i A cp(x-P y*W, y+©e) = ^ Acp (x-0z,y) 


The second integral in equation (l3b) is simplified by introducing 
the notation of the finite part. Define the symbol -f by the equation 

r (x) 5 f (y)dy 1 5a £ a ^ f(y)dy _ 5 f(y)dy 

Jo 5x a/ 1 a-y 2 5x^0 ( a— y) s / 2 5x Jo a/ a— y 

When applied to a single integral this definition is consistent with the 
usual ones for the finite part given as 
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r a f(y)dy = r a f(y)-f(a) f(a) _ _ g ' /*• f « (y)dy + f(o) ' 

Jo ( a— y) 3 / 2 Jo (a-y) 3 / 2 /"T ” U> /iPy + /T - 


However, when applied to double integrals an inconsistency with regard to 
the order of integration between the two symbols | and £ arises. 
Hadamard (reference 2) and Robinson (reference 3) both use the convention 

that the order of integration in the operation j // f(x,y)dy dx is 
reversible; that is, 

/ dy / dx f(x,y) = | / dx / dy f(x,y) 

Such a convention excluded from the area of integration all singularities 
over which the order of integration is not reversible. These singular 
regions are then treated separately. This convention has the disadvantage 
that, in constructing a series of integrals, the value of a given inte- 
gral is not independent of succeeding integrals. 

The operator -f avoids the above difficulty. The value of an inte- 
gral defined by f(y)dy is independent of succeeding operations. 2 At 
the same time, however, the order of integration of operations involving 
the sign -f- cannot be reversed. Hence 

f dy dx f (x,y) ^ / dx -f- dy f(x,y) 


2F or example, according to reference 2 


nt an 
L ( !-O s/ Vn 




However , 

r l d i = r x a5 r l ±1 

(s-n) 0/2 ,/n a 0 (S-n) 3/ ‘VT 
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Return now to the second integral in equation (13b). Applying the 
uefinition already mentioned for j~, it is found for this type of inte- 
gral that 

f(x,T) ay r A4G js_ 

J A(x)-G(x) ^G 2 — (A— y) 2 JM 

By means of the last formula, the equation for <p beeores 
cp = i A<p(xH3z,y)- ~ J dzi J? dyi 

where r c is given by equation (2) . 

Notice that if the integration had been made first with respect to 
xi the result would be 


L ./p2_/fi_TT^2 J 


cLy 



z(x-xi)A cp(x x ,yi) 

[ ( y-yi ) 2 +z 2 ] y(x-xi) 2 -P 2 (y-yx) 2 -P 2 z 2 


which reduces immediately to 


Equations (lta) and (l^b) illustrate the vital importance of the 
order of integration. In fact by subtraction the result can be derived 
that 

/ toi £ 371 - f 471 £ 311 ~ 7 s ^ u ^ z,y) 

If the fj,Tj,z coordinates are used, equation (13) becomes in the 
notation presented in equation (5) 
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The limits of integration depend 
now upon the position of the 
sizes £ and T) with respect to 
the Mach lines in the x,y plane. 
The area t is still hounded "by 
the curve r c 2 =0 and infinity. 

The asymptotes for the curve 
r c 2 =0 are given hy the two equa- 
tions 


Pi , , ms+p 

t»-*U = “ — — r 

(j-2 mi-P 


v4 


Hi t t t ^ “2“^ 


asymptotes The sketch shows how the 

area t in the xy plane trans- 
fers to the |,T| plane. In the 
| mi < P case in which hoth mi and m£ 

m 2 < ]3 are less than P (the case for 

which the sketch was drawn ) , the 

asymptotes are straight lines having positive slopes and the limits of 
integration are always from one side of the cone to the other and from 
minus infinity to some maximum value. Consider the case in which the 
integration is made first with respect to T|i. Then defining 




there results 


z/l d 1 T? 

^ e-£o 2it^n z dq Hi ^1/sLoo 


H 2 v(z^ 


1V) 


2n,l-HBg 




Lq+Li 


L 0 -Lx 


[(i-m) g; + (t-ijjgj * » 
{ [ ( "-' i) s: -<*-Q e ] + z2 } 1 


where the z 2 has been increased by e since again when the 5 derivative is taken through the first 

integral an indeterminate form results s The evaluation of equation ^ 18^ ■nT»nr*na^ a 1nfl+ a a -Pah nnno. 

tion"(l3). Thun'” " " ' '' ' ' ' “ " 


< Ml ) a + 

Hi mi.-HDp 

{["(n-m) A 

L 1 1 Ui m-i+nio z I 


m _ !4. Z 1 P L 0 ,+L l' a _ 

'T ~ •'-■‘-“i / tniX 

0 2jc mJ Lo '-Ia« 


J \ Z H2*/p 2 -^1 2 .. ^ 

( PV^“ Bh-Hno > *11 / 




I^+Lx 


— J. 


(1^) £ + ( S-tJ g_ 


A cp( £x>T|x) 


where the prime on the L indicates its value when £ is given by the upper limit of equation (18). 
The firBt term in equation (19) becomes 

fflxz(mxflfe+ft g ) znig £ P g -dx s 

n™ _z i^x-HUa) *i P g -mx 5 mx+ni2 „ , 3^2 dn-i 


li.m -* y P -m gx^ J A „ / „ „ 

. ^0 2% f f zCmxma^) _ /^g~|a 2 1 ( ^ X 


l L(mx-Wfe)-y p 2 -nix 2 blx-htis J j 


0 1 — JUl _ / f Xoi'Hl l } ^ Tl 1 •— Xi } 


» — i-' * •- 
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which reduces to 


where 


| A cp (g a >T}a) 


/ 2 £ 

_ ZM-g */ P -mi 
£ a “ 5 mi+ni 2 

ziiiCmiuig+p 2 ) 

na = n r . 

(mi+mg) y 3 a -*ii 2 


( 20 ) 


The second term in equation (19) simplifies when the finite— part notation 
is introduced so that finally 


9 = | Aq> (£ a >na) - 


zP 2 (mi+ma) 


f ih £ 


a Acp(g 1 ^ni) 
dTl1 £3— 


(21) 


If the integration had "been made first with respect to the results 

would have "been 


q> 


- 1 A „( s^) - MSai-gi f ini f at 4 ( 22 ) 

2 2rtLLiUg J J T r„ J 


where 


*b - 


n b = n 


zp -1 vP^-ni 2 2 


mi+mg 


ziig(mim2+3 2 ) 


(mi+mg) 


(23) 


The geometric interpretation of the points 5^^% and £ a ,T) a is 
a3 follows: They are the points at which the forecone in the | ,T| 




■P <H U/l 
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There exists the interesting corollary obtained by subtracting equa- 
tion (24a) from (24b) ; namely, that the difference between an integration 
of supersonic doublets made first in one order and then in the reverse is 
equal to the difference in the magnitude of the distribution at two points 
in the plane. 


The Vortex Distribution 


The velocity potential at the point xi,yi*z of an elementary 
horseshoe vortex at the point xi,yi,0 is given by the equation 
cp = — z(x— xi)/2n [(y— yi) 2 +z 2 ]r c . If fl in equation (3) is taken to be 
the induced velocity u in the direction of the free stream, and, as in 
the case involving doublets, the flow field is considered to contain no 
bodies with thickness (so A(du/dz) is everywhere zero), then equar- 
tion (3) reads 


_1_ _6_ P Az(x-xi)Au dxj_dyi 

2n dx J r J. [(y-y i) 2 +z 2 ]r c 


Now, since by definition 



u dx 


the potential for a distribution of vortices can bq written 


(25) 


1 , f P z(x-xi)Au(xi,yi)dxidyx 

2 rt JrJ [(y-yi) 2 +z 2 ] r c 


( 26 ) 


In terms of the £,r|,z coordinates, equation (26) transforms by 
simple substitution to the equality 


cp = 


(mi+xnaJz 

2rtp 1 p 2 


U rn 


<™> £ * <wx) 5 


Audi, Til) 


{[ (tmi) (& "^ } h . + zS } 


dTii 

(27) 
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APPLICATION OF FUNDAMENTAL FORMULAS TO THE PLANE OF THE WING 


Pressure in Terms of Slope for Symmetrical Bodies 


Since in linearized theory the equations for pressure coefficient 
and surface slope are 


C 


P 




> 


J 


(28) 


the equation for the source distribution, equation (ll)> can be rewritten 
as 


C 


P 


2 b_ p jP ^uC ^i , yi ) dxidyi 
« 8zo(-J r c 


(29) 


where is the value of the slope on the upper surface of a symmetri- 

cal body. In terms of the £,ti,z system this becomes (compare the 
transition from equation (13) to equation (15)) 


C =2/^l__3_ + _l_8 N \rr lu( d Ti 

P * \ H2 dq pi WJJ T J r c 


(30) 


In carrying the partial derivative through the integrals there results 
(as in equation (19)) for 0 < mi < (3 the equation 


Cp = lim 


Lq 1 +Li 1 '"u 


-8-r o wa 

e-> 0 ^iOl o ‘-Li« 


In l 5 “ 


zp 2 yp^-mi 2 


mi+ms 




2(mi+m2) 


/ ds i f 4 " 1 




1 (ijjiu) 


and this reduces to the expression 
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for 0 < mi < (3, 0 < m 2 < <» 


C 


P 





2(m.i+xn2) 



<™>;* + (t-tpj* 


r 3 
i c 




(31a) 


Similarly, for 0 < m 2 < P, 0 < mi < » 


Cp Tp^te* 






H2 




(31b) 


and for f3<ma <°°, P< m 2 < » 


C 


P 


/ \ m i / \ ■ m 2 

2(mi+ma) r jir + ilg 

rt|aifi 2 J J T Tc a 


Mi^i) dg x dt]x 


(31c) 


where Sa'^a'^b' an<3 ' % are given by equations (20) and (23) 

As special cases of the results given by equations (31), consider £ 
and q to represent the x,y and the r,s axes, respectively. In the 
former case, x replaces £ as ms — 00 and y replaces q as 
mi 0; hence 


c p = jr Mx-Pz,y) - f^dxi jfdyi — — > ~ u - 3 Xx,yi) ' (32a) 

and — 

C P = -| /Vi Xto, fesdishiuii (32b ) 

J J T r c 

In the latter case r replaces £ as ma-> p and s replaces q as 
mi — > P; hence 
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C P /y (£Sil±^2lI x u (rx, Sl ) axxd Sl (33) 

T r c 

In each, of the equations (3l) through (33)# z can he set equal to 
zero without affect ing the validity of the equalities. When z vanishes, 
the following identities hold 



With these equations the pressure coefficient on the surface of a sym- 
metrical n onl ifting wing can he determined if the surface slope is given. 
The special cases of equations (3l) in the plane z=0 and in the x,y 
and r,s coordinate systems are given in the summary of this section 
(see equations (4-7) and (50)). 

Vertical Induced Velocity in Terms of the Jump in Potential 

It is proposed next to find the vertical, induced velocity in the 
z=0 pi ane as a function of the jump in potential across that plane. 

Consider equations (24) for the doublet distribution and take the 
partial derivative with respect to z of both sides; then find the limit 
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of the resulting expression as z goes to zero. If equation (24a) is 
used, for example, there results for the first term 



which becomes 


2(mi+ma) J 0 2 -mi 2 4 
and for the second term 


M-s(P S — mi 2 ) Acp(| ,n)-Hil(mim2-H3 2 ) ~Acp(£,Ti) 


(36) 


-0 


A<p(li,m) 


r C 3 


But this reduces to 


^(mi+ma) 

2nni|i 2 


'faff* 


A Cp(|i,Tu) 


P g (m!+ma) _i±% p 
23t|i 1 |i 2 



A<p(Si,tji) 


(37) 


and since 



dm 

[(Xg-T)i) (t^i-Ai )] 3 / 2 


= 0 


the second term in expression (37) vanishes* -Finally, therefore, the 
vertical induced velocity in the plane of the wing w u becomes 
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v = — 
u 


2(1311+02) y 0 2 -mi 2 




Acp (|^)+Hi(mim2+P 2 )^- A?( £,ti) 


P 2 (mi+m2) 

2khi|i 2 



Acp(£i,ru) 


r 3 

r o 


and similarly 

for 0 < ms < P, 0 < mi < 00 


( 38 a) 


w u .= - 


2(mi+m 2 )*/p 2 -^2 2 L 


Mi(P 2 -hds 2 )^- Acp(|,Ti)+(i2(mim2+p 2 ) : ^- Acp (|,ri) 

8t} 


0 2 (mi+mg) 

2«nnx a 



«1 


Acp ( |i>t}i) 


Y* ® 

r o 


(38b) 


and for 0 <102 < 0 < mi < 00 


0 2 (mi+m2) rr 

213 ^. 1 ( 1.2 J J . t 


Aq>(|i#nx) 


d£i dT]i 


(38c) 


The special cases of equations (38) obtained when the |,r; coordi- 
nates represent the x,j or the r,s systems are given in the summary 
of this section (see equations (48) and (5l))« 


Vertical Induced Velocity in Terms of the Loading 


The equation for the loading coefficient in linearized theory can 
be written 


Ap 2Au 

T = ~ 


(39) 
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so that equation (27) for the vortices yields for w u the expression 


w 

u 


Ijto 8 (mq.+ma)z7 0 
z 5? 4rt|iiP2 



(tmu) sa 

Hi 1 H-2 




^ (l!^l) 


(Vai) ~ - U-4i) -p 

Hi H 2 


+ Z‘ 


d| drii (4o) 


The evaluation of w u can be divided into two steps; first, the pro- 
cedure necessary for carrying the derivative through the first integral, 
and second, the calculation of I where 


I 


lim 
~ z->0 


(mi+nfe)y o w 

' 1 n ' ” X 

W1H2 







+ z' 



dfca. dru(la) 


Again the order of integration is important. To begin with, the 
first integration will be taken with respect to qi. Further the case 
0 < mi< (3 will be considered. Hence, the equation for w u becomes 
(just as in the derivation of equation (l8)) 


w u - z — >0 


(mi+mg)V 0 8 


e — 5>-0 . 3z 


/; 


P 2 '/z 2 +e 2 Vp 2 — Ell. 2 
mi ■fins - 


pho+hx 

J L 0 -Li 


( Ml ) 2i + (R) - 

' 1 Hi ' b 1 Uo 


H2 J 






(ti-tii) — - (i-i x ) — 


H-2 


+ Z‘ 
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where the values of L 0 and Li are given by equation (17) and the 
limit as e goes to zero is to be taken first. Equation ( 36 ) reduces 
to 


llm T 0 ,V f3 2 -mi 2 ^o'+Li 

; ->0 , / 

-^•0 4it|ii \/z 2 +£ 2 ^L 0 '— Li 


(■n— ^11) ^ + 




p p 

zms */ (3 ■ -mi 
itti+ms 


(ia.m) 


+ I 


(43) 


("H — Til) 

Pi mi+ms 


+ z-= 


The quantity within the integral of the first term of equation (43) is 
the same as the similar term in equation ( 19 ) ; hence by analogy 


v. 


u 




+ I 


(44) 


The evaluation of I requires some care. Consider the following 
integral which contains all the difficulties involved in I. 


lim 


/A 

zf(yi,z) 

Ja dz 

-z 2 +(y— yi) 2 _ 




b > y > a 


where f(y, z) and its derivative is bounded and continuous in the inter- 
val a < y < b. 

Integrating by parts 


T - - 3 - im .A. 

o " z dz 


y— 0, y— ^ 

f(a,z) arc tan f(b,z) arc tan + 

z z 


nb 


df(yi>z) ___ y-yi. 

— *. arc tan dyi 

oyi z 
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and since 


lim 
z ^ 0 


dfz(yijz) 

dyi 


arc tan 


y-yi 

z 


dyi 


«fz(y#°)- 


it 

2 


f z (a,0)+f z (b,0) 


this becomes 


I 


o 


f(b) _ f(a) 
y-b y-a 



y-yi z 


It is, however, more convenient to write I 0 in terms of f(yi) 
and not its derivative. In order to do this, it is necessary to intro- 
duce another notation involving singular integrals. The concept of 
Cauchy 1 s principal part is adopted and defined in the following way: 



f(yi)<3yi 

y-yi 



y-yi| dyi 


This procedure can be generalized (see also appendix D in reference 4) 
so that for 



f(yi)ayi _ P f(yi)ayi 
(yi-y ) 2 dy y x -y 


d_ r h f(yi)dy x 
dy '-'a yi-y 


(45a) 


This definition of the symbol f can be made in another way. If the 
Indefinite integral is expressible in the form 


then 



£(yi)<3yi 

(yi-y ) 2 


= G(yi,y) + constant 


— ^ = G(bjy) - G(a,y) 

(yi-y) 


(45b) 


Thus the conventional rules for evaluating a definite integral from the 
indefinite integral can be used. 

Using equation (45b), one can derive by integration by parts the 
relation 
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f* f(yuO)dyi _ f(b‘) _ f(a) _ P b f « (yi,0)dyi 

J a • (yi-y ) 2 r y-* y-a J a 


Hence the equation for I n becomes 


T f(yi»0)dyi ^ ^ 

io = 4 ^ — + « fz (y>o) 

^a (yr-y) 


By means of the concepts introduced by equations (V?), I can be 
evaluated 3 and f in ally v u can be written 

for 0 < mi < 0, 0 < nfe < 00 


a _ _ a ( + 


4 i 


nii+nte 


/ a 5l / a ni 


(^ 1 ) + (S-gj.) g- Ap/q 


Similarly, for 0 < < 0> 0<mi<eo 


T— ■— (§>H) + 


nii+njg 

WllJ-2 


/<*/ ^ 


and, f inally, for 0 < mi <oo > p<m 2 <co 


V-^ mi+ni2 

V n WlU2 


// 


(1-U) £ - (l-l,) £ 4P/1 


' £ 


- (&-i x ) 


1 l 2 r. 


d|i dm 


3 In the evaluation of I the term representing f(yi>z) can be written 
as a function of z 2 . Hence, the term rtf z (y,0) in equation (^-5c) 

vanishes. 
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Again the special cases of equations .(46) when g,T} "become x>y 
or r,s are given /in the following summary X see equations (49) and 
(52)). 


Summary of Formula for w and Cp in the Plane of the Wing 


For the x ? y coordinates .— In these cases g-^x, y m.x-^0, 

and m 2 00 . 

n - a U W ) - a /W an (1.7a) 

P P u «j ■ J T [ ( x— xi ) 2 -P 2 ( y— y 1 ) 2 ] s ^ 2 


C P 


(x-xi)\ u (xi,yi) 


T [ (x-Xx) 2 -p 2 (y-yi) 2 ] 


. 2 , 3/2 


(47b) 


w u = - | Au(x,y) -^-Jdxx jf ayi 


A(p (xi,yi) 


[ (x-xi) a -P 2 (y-yi) 2 ] 3/2 


(47c) 




A t? (xi,yi) 


[( x-xi ) 2 -P 2 ( y-yx ) 2 ] 3/2 


(48) 


£ - - 1 h / to1 / 371 


(x-Xx) (Ap/q) 


T ( y-yi )V ( x-xx ) 2 -P 2 ( y-yx ) ‘ 


(49a) 


vu 1 r, r , (x-xx 


(x-xx) (Ap/q) 


(49b) 


*x) -p (y-yx)‘ 


For the r,s coordinates .— . In these cases g -^-r, r| ->s, mi -=*0, 
and m 2 — . 

C rf (^xHU^si) \ 11 (rx,s 1 )dr 1 ds 1 

p 2*P JJ [(r-rx) ( s— si) ] 3 / 2 


( 50 ) 
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Mo r f&y (rijSi)dridsi~ 

u 8rtJ , J [(r-n)(s— si) ] S/ ^ 2 

% _ r r [ (p~^e*i) t( &— 3i ) 3 (Ap/q,)dridsi 
v o [( s -si)-(r-ri)] 2 y (r-ri)(s-si) 


(51) 

(52) 


It should he emphasized that in the r,s coordinate system the order of 
integration is immaterial. 


S0M1 ALTERNATIVE EXPRESSIONS 


It is somet ime s very convenient to be able to eatress the equations 
given in the previous sections in a slightly different form. Consider, 
for example, equation (ll) which gives the velocity potential due to a 
distribution of sources, thus 


1 P rw u (xi,yi)dx;L<3yx 

77 

where the area t can be defined 
( see sketch) as the area bounded 
by the wing plan form 4 and the 
forecone from the point P(x,y, z). 
Define now the area as being 

the area bounded by the plane 
x=xi and the wing plan form 
ahead of the plane x=xi ( see 
sketch) . It is apparent that r c 
is a pure real quantity everywhere 
inside the area t and is a pure 
imaginary qua n tity everywhere out- 
side r and inside t 1 . The same 
is true, of course, of r c 3 . It 
is clear that all other terms in 
the integrands of the integrals 
which have been considered in the 
preceding section, in particular 
■WuCxijTi) in equation (ll), are 



|Xi 

4 The actual definition of ~ that it Is the area within the Mach fore- 
cone, is often replaced by the one Used here since the strengths of the 
sources, etc., are zero ahead of the wing. 
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always real over the entire wing plan forms, Hence the integral 


1 P P v 

*Jr j 


is a complex quantity, the real part of which is the velocity potential; 
thus 




1 P p w u (xi,yi)axidyi 

it J J r c 

T i 


Similarly, each of the integrals in equations (47) and (48) may he 
replaced by the real parts of their values taken over the area t x . 


(53) 


The evaluation of the terms involving the finite part are particu- 
larly simple when the area is used since, if 


1 f(y)dy 
(*-y) 3/2 


= F(y)+C 


then for positive a and b ^ a 

r b f(y)dy 

r.p.-f 


Jo ( a-y ) 3 / 2 ' 


= r.p. [F(b) -F(o)] (54) 


For example , consider 


i 7237 


J o (aa.ysjs/a 


r b ^ay 

r.p. 4- - — 

(a 2 — y 2 ) 372 


where a < b. From the relation 


/ 


y g ay 


— arc sin 


(a 2 — y 2 ) S//S */ a 2 — y 2 


together with equation (54) 

I = r.p. - y- £-t- ; - 

\ y a 2 -b 2 


• b \ it 

arc sin — = — rr 

a. J 2 
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Notice the simple extension that 

1 7 2 3y p'b ^dy 

= r.p. -L = — jt 

(a 2 — y 2 ) 3/2 ih (a 2 -^) 3 / 2 

Further examples using the real part of the integration taken over 
the area will he given in the next section on special applications. 



SPECIAL APPLICATIONS 


An Integral Equation 


Applications of the results given by equations (48) through (52) are 
apparent. One of the more important uses, however, comes in the develop- 
ment of integral equations necessary for the solution of many supersonic 
wing theory problems. 


An example of such an application 
arises in the analysis of the slender 
rectangular wing at an angle of attack 
a. Since the wing chord is long coxa- 
pared to its span, and since along the 
side edges the loading falls to zero, 
an approximation to Ap/q is given by 
the equation 

££ = 4af(|) yi-<y/s) 2 (55) ■ 

where f(x/s) is an unknown func- 
tion. The assumption made when 
using equation (55) is, of course, 
to fix the spanwise variation of 
loading but leave the chordwise 
variation arbitrary. 



The function f can be determined by the condition that the value 
of w u in equation (4-9a) is a constant all along the center line of the 
wing. The area r is indicated. by the shaded region in the last sketch 
so that equation (4-9a) for the case y=0 (and for added simplicity 
^=l) becomes for x > s 
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^-s 7i 2 ,/ (x— Xi) 2 — yx 2 


I T* (z-n)f('n'W l! f I 11 7 i -(yi/e) a _ayi 

A-b V s J yiV (x - Zl ) £ - yi 2 


and for 0 < x < s 


*—(!)♦ I/* W (^>/ 

V ° o V / Mxhxi) 7i v (2— Xi ) 2 — yi 2 


Introduce the notation 


Sn. = TT> 9 - f . ki = ka = 8 - ^ 


and these equations become , 5 since a = — 

v o 


for 0 < 6 < 1 


for 1 < 9 


1 = f(0) + |y 9 kaB 2 f(0i)dex 


O p ® O P 

f(e) + | / k^B a f(e 1 )de 1 + | / ExfC ei)aei 

Jfl_n Jr, 


3 The symbols B and E indicate elliptic integrals. Thus 


E n =/ 


'3- /l-k 2 t 


■k n 2 t 2 ' p 1 dt 

dt; % - / ; ■ ; «M 

Jo J (l-fc a )(l-tg t=) 


Bn = 


E n -U-k n 2 )Kn 
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The solution to these equations has been obtained numerically ann 
has the form shown in the sketch. 


Since from equation (55) the 
average spanvise loading P can 
be readily calculated as 

l) 

the curve for average span loading 
can easily be constructed. This 
curve is also shown in a sketch 
together with a portion of the 
variation of P obtained from an 
exact linearized analysis. 

In the interval where the compar- 
ison can be made (i.e., near the 
leading edge) the agreement will 
be the poorest because in this 
region the spanwise variation 
deviates most radically from the 
value as siimed in the construction 
of the integral equation. 


f(e) 




Drag Reversibility Theorem 


The well-known theorem that the drag of a symmetrical nonlift ing 
body is the same in forward and reversed flight at the same speed (see 
reference 5 or 6 ) can be derived in another way using the results of the 
preceding sections. 

By definition 

CD = S Jg f 2 X u(x,y)C p (x,y)dy dx ' ( 57 ) 

where 

S • area of the wing 
X u slope of the upper surface 
Cp pressure coefficient 

Using equation (^7b)j 
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dx 



dxi 


(x-xi)X u (xi,yi) 


( 58 ) 


How the equation for the drag coefficient in reversed flow can he 
obtained by: 

1. Replacing the area by t 2 such that r x + t 2 = S 

2. Rotating the axial system in the xy plane through l80° 

3. Reversing the signs of X u (x,y) and X, u (xi,yi) 

There re stilts 


— / / Mi(x,y)dy dx 


r.p 




dxi 


(xi-x)X u (xx,yi) "I 


( 59 ) 


and subtracting equation (59) from (58) gives 


Cp-C Dr = “ ~ r.p.y’ dyy° dx J' dyijT dxi 


gWxijyi)M*jy) 




xi*. u Ui,yi)M x >y) 


37 ® 
r C 


(60) 


Since the symbols ±i>yi#x,y are dummy variables of integration, the 
last term in equation (60) can be written 

sr r - p -/ s ^7 — 

and reversing the operators / dyi /dxi and / dy / dx (but always 
preserving the same order within the operation) yields for the second 
term in equation (60) the same expression as the first term except for 
sign. Hence 


_ °D r = 0 


or 

°D = °D r 


as was to be shown 
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Lift on Wings With. Supersonic Edges 


The lift on any wing can "be written 


L 

q 



dy ds 


Moreover, 





where T.E. and L.E. denote 

the trailing edge and leading edge, 

respectively, and t.E 

value of the velocity potential on 

the upper surface of the wing at 

the trailing edge. 

Consider now a wing with al 1 
edges supersonic and a straight 
trailing edge not necessarily at 
right angles to the free— stream 
direction. Let the wing he a 
plate having arbitrary twist and 
camber. Then for, a point on the 
wing, the velocity potential from 
equation (53) can be written 


t 



9 




v u(^x/yi)dxidyi 
( x— x 1 ) £ — {3 2 ( y— y 1 ) 2 


and if the equation of the trailing edge is 


x = a + y tan A 

where a is some constant, then 


'PT.E. 


ir.p. 


v u (xi,yi) dxi dyi 


S ^ J (a+y tan A— Xi) 2 -P 2 (y— yi) 2 


so that the total lift on the wing can be written 
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L 

i 


k 

*v 0 


r.p. 



w u (zi,7a)dJtidyi 


a +y tanA-zi ) 2 -£ 2 (y-yi) 2 


The area S, being that of the wing plan form, does not depend on y 
so the y integration can be made first and, since the edges of the wing 
are supersonic, the interval si < y < S2 must always contain the roots 
Xi and \ 2 of the expression under the radical. Hence 


r.p 


■L 


s 

— s 




rX 


X 2 


<3y 


J (p 2 - tan 2 A)(X, 1 -y)(y-A 2 ) *a J (p 2 - tan 2 A )(\i-y) (y-^. 2 ) 


and since 


then 


■r 


dy 


= it 


L 

q 


-4 


7 (>-i-y)( 7-^,2) 

r r w u( x i,yi) 

JsJ v 0 


J P 2 — tan 2 A 

An alternative expression for equation (6l) is 

43 


dxidyi 


(61) 


C L " 


y P 2 — tan 2 A 


(62) 


where ct, is the average angle of attach of the surface and by defini- 
tion 

--ip r w uUi>yi) „ , 


It is interesting to notice that the lift coefficient for such a wing is 
the same__as that for a two— dimensional flat plate flying at an angle of 
attack a into a free stream, the speed of which is given by the com- 
ponent of velocity normal to the trailing edge of the three-dimensional 
wing just studied. This result has been derived previously in refer- 
ence 7. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Oct. 1 6, 1950. 
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